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ABSTRACT 

The  detection  of  the  Zeeman  effect  by  spin  precession 
experiments  has  led  to  the  development  of  magnetometers 
which  utilize  vapor  or  gaseous  sensing  devices  capable  of 
high  precision.   A  recent  example  of  this  type  of  instru- 
ment is  the  helium  gas  magnetometer.   The  purpose  of  this 
paper  is  to  present  the  background  theory  of  a  helium  mag- 
netometer which  utilizes  a  modulated  pumping  light,  and  to 
describe  the  equipment  used  to  construct  an  experimental 
magnetometer  system.   The  feasibility  of  this  system  as 
an  instrument  to  measure  changes  in  the  earth* s  magnetic 
field  is  confirmed  by  comparison  with  a  commercial  magne- 
tometer. 

The  author's  work  on  this  device  was  accomplished 
at  the  Varian  Associates  research  laboratory  at  Palo  Alto, 
California,  during  the  period  January  to  March,  1961. 

The  author  wishes  to  express  his  appreciation  for  the 
assistance  and  encouragement  given  him  by  Drs .  A.  L.  Bloom 
and  W.  E.  Bell  of  Varian  Associates  during  this  period. 
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I .   Introduction 

The  utilization  of  the  earth's  magnetic  field  for  man's  needs 
can  be  traced  back  through  many  centuries.   To  further  his  advantage, 
man  has  undertaken  study  of  this  phenomena,  and  it  remains  one  of 
today's  challenging  geophysical  frontiers.   The  present-day  require- 
ments for  sensitive  underwater  detection  devices  have  increased  the 
need  for  precise  magnetic  field  measurement.   To  meet  these  needs 
for  measurement,  a  class  of  instruments  has  been  developed  which 
uses  vapor  or  gas  as  the  detecting  substance.   The  underlying  princi- 
ple of  their  operation  can  be  summarized  as:  (a)  accurate  detection 
of  the  Zeeman  effect  within  the  working  medium  and  (b)  subsequent 
utilization  of  this  inforrnation  to  form  a  self-correcting  closed 
loop  system.   In  this  manner  the  exterior  magnetic  field  information 
may  be  extracted  when  desired  or  on  a  continuous  basis  by  suitable 
recording  devices.   Such  magnetometers,  utilizing  rubidium  vapor, 
are  commercially  available. 

The  work  with  helium  gas  as  a  magnetic  field  detection  device 

has  been  largely  directed  towards  obtaining  usable  resonance  signals. 

1  2 
However,  operating  magnetometers  using  helium  hav-  been  reported. 

The  principle  of  optical  pumping,  together  with  resonant  mixing  of 

the  metastable  helium  sublevels,  provides  a  means  of  detecting  the 

magnitude  of  the  Zeeman  effect.   This  action  can  be  visualized  as 

one  which  allows  the  pumping  light  to  cause  an  unbalance  in  the 

metastable  population  densities  followed  by  a  redistribution  under 

the  control  of  external  electronic  devices.   Fortunately,  this  cycle 

P.  A.  Franken  and  F.  D.  Colgrove,  The  Metastable  Helium  Magnetometer 

(1) 

2 
Magnetometer  System  for  Orientation  in  Space,  Electronics, 

April  8,  1960,  pp.  55 


can  be  detected  easily  by  monitoring  the  residual  pumping  light,  i.e., 
the  output  light  which  passes  through  the  sample.   This  output  light 
will  contain  information  concerning  the  process.   Both  polarized  and 
unpolarized  light  are  used  for  this  purpose.   The  unpolarized  light 
causes  a  different  distribution  in  the  population  densities  from  that 
caused  by  the  polarized  light.   In  the  case  of  polarized  light,  a  net 
magnetic  moment  is  formed  by  increasing  the  density  of  one  sublevel 
at  the  expense  of  the  other  sublevels  involved.   As  a  practical  matter 

this  difference  leads  to  dissimilar  physical  effects  and  no  detailed 

3 
comparison  will  be  made  here. 

The  desired  information  contained  in  the  output  light  can  take 
several  forms  depending  on  the  method  of  operation  desired.   In  the 
simple  detection  scheme  described  above,  the  f£^itput  light  information 
will  be  in  the  form  of  a  resonance-like  amplitude  variation  which  has 
a  peak  position  that  corresponds  to  the  strength  of  the  ambient  mag- 
netic field.   This  variation  can  be  achieved  by  1)  introducing  into 
the  gas  a  radio  frequency  field  which  corresponds  to  the  Larmor  pre- 
cession frequency  and  then  2)  varying  the  radio  frequency  or  the  total- 
magnetic  field  in  a  periodic  manner.   As  the  changing  frequency  or 
magnetic  field  passes  through  resonance,  a  change  in  the  light 
absorption  will  take  place. 

A  second  method  of  utilizing  the  pumping  light  is  found  by  fix- 
ing the  applied  radio  frequency  field  at  the  Larmor  precession  rate 
and  examining  the  radio  frequency  content  of  the  output  light.   It 
will  be  found  that  the  output  light  is  modulated  at  the  applied  radio 
frequency  and  a  closed  loop  oscillator  can  be  formed  at  the  Larmor 

•a 
For  the  unpolarized  light  physics,  see  F.  D.  Coi grove  and 

P.  A.  Franken  (2) 


precession  frequency.   Alternately,  it  is  possible  to  heterodyne  two 
radio  frequencies  which  are  within  the  line  width  of  the  Larmor  pre- 
cession frequency  and  extract  a  difference  frequency  from  the  light 
for  magnetometer  purposes. 

The  magnetometer  constructed  in  the  course  of  this  investigation 
utilized  a  principle  of  "optical  modulation"  or  as  alternately  phrased*, 
"pumping  light  modulation".    This  method  utilizes  a  polarized  pumping 
light  to  change  the  population  densities  of  the  metastable  levels  and 
a  pumping  light  modulation  at  the  appropriate  Larmor  precession  fre- 
quency.  Passing  the  light  through  a  helium  sample  results  in  a  net 
spin  orientation  within  the  gas.   Detection  of  this  orientation  is 
accomplished  by  monitoring  the  pumping  light  intensity  after  it  passes 
through  the  helium  sample.   The  amount  of  absorption  will  be  altered 
as  the  modulation  frequency  passes  through  the  natural  spin  precession 
frequency  corresponding  to  the  ambient  magnetic  field.   This  latter 
method  was  chosen  in  preference  to  the  previously  mentioned  methods 
because  the  theoretical  conception  is  of  recent  origin  and  not  necess- 
arily because  it  is  the  most  desirable  as  a  practical  system.   It  was 
the  intent  of  this  investigation  therefore,  to  show  the  feasibility  of 
a  magnetometer  using  modulated  light  and  compare  its  performance  with 
an  existing  system,  in  contrast  to  making  a  detailed  comparison  with 
other  possible  helium  systems. 

In  order  to  present  the  results  of  these  efforts  in  a  comprehensive 
manner,  this  thesis  contains  a  brief  discussion  of  the  theoretical  back- 
ground as  aided  by  subsidiary  proofs  contained  in  the  appendices.   This 
is  followed  by  a  discussion  of  the  theoretical  magnetometer  design 
which  subsequently  materialized  into  experimental  hardware. 

A.  L.  Bloom  &  W.E.  Bell,  Radiofrequency  Mixing  in  Optical  Pumping 
Experiments,  to  be  published. 

SW.  E.  Bell  &  A.L.  Bloom  (3) 


II.   Theoretical  Background 

A.   The  Helium  Energy  Levels 
I.  The  Fie Id- free  Spectrum 

From  the  viewpoint  of  classical  mechanics,  the  atomic  electron 
motion  can  be  analyzed  in  terms  of  certain  assigned  parameters.   The 
first  of  these  arises  from  the  motion  of  the  electron  in  relation  to 
the  nucleus  and  is  expressed  in  terms  of  the  orbital  angular  momentum 


which  has  a  magnitude   A  JP(0  +  1)        :r  ^=       where  J(      is  the 
azimuthal  quantum  number  and  h  is  Planck's  constant.   To  a  first 
approximation,  this  parameter  can  be  applied  to  a  complicated  atom 
by  assigning  to  each  electron  an  individual  orbital  angular  momentum 
and  taking  their  vector  sum  to  represent  the  total  orbital  angular 


momentum  JlJ^/^+j\       ~   A  t-  Since  analysis  by  quantum 

mechanics  will  show  that  the  individual  orbital  angular  momentum  can 
only  be  oriented  in  certain  discrete  directions  to  one  another,  the 
vector  sum  can  only  assume  certain  values,  namely,  L  =  0,1,2,.... 

In  a  similar  manner,  the  assumed  rotation  of  the  electron 
around  its  own  axis  will  result  in  an  electron  spin  angular  momentum, s 
This  combines  with  other  electron  spins  to  give  a  resultant 
-£-  i/</'(j.A  =  Ajr.     Again,  only  certain  discrete  values  are' 

3IM/ $(:>  +  /;     ^  .   "* 

hi       H      i 
allowed-  S  taking  on  the  magnitudes   3">~5"'»;  -  -  -  -  Q 

where  N  is  the  total  number  of  electrons. 

The  total  angular  momentum  can  be  found  by  the  vector  sum  of 


the  orbital  and  spin  momenta  to  form  JL-  \   J~(%T*^  u     ~™   — ~~~ 
where  again,  the  discrete  magnitudes  and  directions  result  in  the 
rule:  "The  vector  sum  of  S  and  L  is  such  that  the  different  possible 
values  have  integral  differences".   It  shall  be  considered  here  that; 

J'  =  L  +   S 

A 

defines   the   total   angular  momentum  vector  .in   terms  ot       rjfrp    » 


The  helium  spectrum  is  composed  of  singlet  and  triplet  lines. 
This  duality  arises  from  the  two  electron  spins,  since  each  has  a 
quantum  number  s  =  %  and  they  can  be  aligned  parallel  or  antiparallel, 
As  a  result,  the  spin  angular  momentum  can  be  zero  or  one: 


s. 


JL. 
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SA 


Fig.  1.   The  Helium  Spin  Alignment 
For  the  S  =  0  combination  of  electron  spins,  only  singlet  lines 
can  appear,  and  the  result  is  designated  parahelium.   For  example,  the 
P  (L  =  1)  state  results  in  a  total  angular  momentum  J  =  1  since  S  =  0 
and  only  one  P  state  per  n  quantum  number  will  appear. 

Orthohelium  is  the  result  when  S  =  1,  for  then  additional  energy 
levels  can  appear.   For  the  P  state  (L  =  1)  three  lines  will  result. 
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Fig.  2.  The  Orthohelium  Angular  Momentum 
Also  for  the  orthohelium  S  state  (L  =  0) : 


s  =  t=  I 


Fig.  3.   The  Helium  S  State  Angular  Momentum 


The  general  quantum  selection  rules  for  transitions  from  one 
state  to  another  are: 

a.  At   =  ±1 

b.  AL   =  ^*>° 

c  AS  -  ^ 

d       A  J   nil     O    [W£     T-o     Cc    T~c~l 

e.  £M  =  ±/,0  £^  M-^  r°  M'°     *-*   &JsC\ 
Rules  a.  and  c.  prohibit  normal  transitions  from  the  ortho  to 
parahelium  states,  and  hence,  such  a  transition  must  arise  from 

collision  (thermal)  excitation.   This  fact  creates  a  metastable 

3 
condition  at  the  S  level  which  is  dependent  on  thermal  excitation 

to  decay  to  the   S  ground  state. 

The  portion  of  the  helium  spectrum  which  is  pertinent  to  the 

magnetometer  application  is  centered  about  this  metastable   S  state 

which  is  uniquely  characterized  by  permissible  photon  excitation  to 

higher  energy  levels  and  thermal  relaxation  to  the  ground  state. 

In  diagrammatic  form: 

% 

3fi 


a  — 7T- 
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Fig.  4.   Helium  Spectrum  in  the  Absence  of  Magnetic  Field 


2.   The  Zeeman  Splitting 

The  rotation  of  the  electric  charges  within  the  atom  leads 
to  the  production  of  a  magnetic  moment  which  is  proportional  to  the 
angular  momentum.   The  effect  of  the  magnetic  moment  is  to  cause  a 
precessional  motion  when  the  atom  is  immersed  in  a  magnetic  field, 
and  in  accordance  with  quantum  mechanics,  only  discrete  directions 
in  relation  to  the  field  are  allowed.   This  means  that  a  "space 
quantization"  of  the  total  angular  momentum,  J,  exists  in  a  magnetic 
field  environment,  and  that  the  total  angular  momentum  in  the  direction 
of  the  field  will  take  on  only  discrete  values.   It  can  be  shown  that 
there  are  2J+1  such  values,  hence,  there  are  2J+1  energy  sublevels 
for  each  field- free  state.   The  actual  energy  associated  with  each 
sublevel  depends  on  the  product  of  the  magnetic  field  intensity  and 
the  component  of  the  magnetic  moment  in  the  field  direction. 

In  the  case  of  orthohelium,  the  so-called  anomalous  Zeeman 
splitting  is  observed.   Because  of  the  Lande  g- factor  involved  with 
this  effect,  the  mathematical  consideration  of  the  Zeeman  energy  is 
developed  in  Appendix  A.   In  diagrammatic  form  the  Zeeman  splitting 
can  be  shown  as  in  Figure  5. 

N.  F.  Mott,  Elements  of  Wave  Mechanics,  Cambridge  Press, 
1952,  pp.  79,  Herzberg,  (M-)  pp  106. 
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Fig.  5.   Helium  Energy  Levels  Showing  Zeeman  Splitting 


3 .   The  Transitions 

As  shown  in  Appendix  C,  there  are  18  allowed  transitions  be- 
tween the  -*S  and  ^P  states.   The  wavelengths  of  these  lines  are  in 
the  near  infra-red.  (1.08  microns)   In  the  presence  of  external 
photon  excitation,  there  are  essentially  four  separate  actions 
which  effect  the  transitions  between  these  two  states.   These  are: 

(a)  Induced  absorption  in  the  presence  of  the 
pumping  light  which  cause  a  transition  from  one 
of  the  ^S  sublevels  to  one  of  the  3p  sublevels . 

(b)  Induced  emission  in  the  presence  of  the  pumping 
light  which  cause  a  transition  in  the  opposite 
direction  to  (a)  . 

(c)  Spontaneous  emission  from  the   P  levels  to  the 
3S  levels. 

(d)  Thermal  relaxation  effects  which  effectively 
retire  and  replace  metastable  atoms  in  all 
levels . 

Of  these 9  the  second  is  generally  considered  to  have  little 
influence  on  the  problem  due  to  the  short  lifetime  of  the  atom  in 
the   P  state „  (a/  10~°  sec.)   Each  of  the  remaining  effects  have  a 
finite  individual  probability  of  occurrence  for  each  level,  although 
the  thermal  effects  are  usually  considered  equally  probable  for  all 
3S  sublevels.   The  lifetimes  of  the  metastable  atoms  have  been 
measured  and  fall  in  the  range  of  10"^  to  20x  10 "^  sec. 


As  shown  in  Appendix  B,  the  effect  of  using  a  circular  polarizer 
is  to  limit  the  allowed  transitions  to  a  total  of  six3  which  are 
those  involving  a  J^  H  of  plus  or  minus  one  depending  on  the 
direction  of  polarization.   This  in  turn  radically  changes  the  sub- 
level  population  densities  as  shown  in  Appendix  Cs  and  as  a  further 
result,  a  net  magnetic  moment  is  formed.   It  is  this  net  atomic  mag- 
netic moment  that  will  allow  analysis  of  the  system  ensemble  on  the 
basis  of  a  magnetization  vector  per  unit  volume  which  will  follow 
under  the  consideration  of  signal  development. 
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M-.   The  Larmor  Precession  Process 

In  Appendix  A  the  Zeeman  sublevel  difference  is 
considered  from  the  viewpoint  of  an  energy  relation.   It  is  of 
equal  importance  to  consider  the  precessional  effects  associated 
with  this  energy,  for  each  is  related  to  the  other  by  the  familiar 
relation: 

AE  -  h  * 


(i) 


From  Appendix  A,  equation  (10): 

where   yet,    -  atomic  magnetic  moment 
-  Bohr  magneton 
<j     -  Lande  g- factor 
i'(-Thi')     —  net  atomic  angular  mementum  in  units  of  ™ 

For  convenience  let  5"-  Qj  /j"(iJ"*  /)    be  called  the  "spin"  vector. 

Then:    ^  s  _^§- 

But:  ^ 


(J it)  — r~   —  °  y    the  gyromagnetic  ratio 


Hence: 


>*  =  -y<  *  "  JF 


And  since: 

Torque  =    ^^    X  /"^ 


7     = 


/  *  ^5  A  // j 
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But,  torque  =  J    (angular  momentum) 

so:   T  =  .L.(ts) 

and:    iE  -  -y(s*h) 

which  can  be  written  as  a  scaler: 

Let    the    spin      j      rotate    around  /•/       an   increment  d5    corresponding 
to   an  angle  <a  y    .      Then: 

a  (3    ^        -f*S - 


and: 

^  =  angular  velocity  -  <IL    =  ^JJi±^L^^     «  Y',-/ 

Hence,  for  a  magnetic  moment  of  ^/^    ,  there  is  a 

precession  around  the  H  field  at  an  angular  rate  £^  . 

If  the  field  is  //,    then  there  is  a  unique  rate,  CUC 
associated  with  the  precession. 

But,  note  that  this  frequency  corresponds  to  the  energy  of 

the  precession  as  developed  in  Appendix  A,  where  &  M   ~  I  , 
Equation  (12)  there,  was: 

For  the   S  levels,  a  short  calculation  will  reveal  g  =  2, 
and  from  (1):     J^    -    Ajtjj    (AiA)      =  Q.jAd 

#  £623  773  X /€>-*? 

where  1  gamma  =  10  gauss  . 
Since  the  earth's  field  is  nominally  one-half  gauss,  the  resonant 
frequency  used  in  the  helium  magnetometer  is  1.4-  mc „ 

12 


B.  The  Signal  Development 
1.   The  Pumping  Process 

It  is  of  interest  to  consider  a  rotating  frame  of  reference 
Let  it  be  assumed  that  the  spin  vector  is  standing  still  in  a  fixed 
frame  of  reference,  and  the  observer  is  rotating  at  angular  rate  UJ 
To  the  observer  the  spin  vector  will  appear  to  be  rotating  backwards, 
and  it  will  move  through  the  angle  A  f?   in  the  time  ZT    .   From 
Figure  6,  it  can  be  seen  that  ZSS  »  &0  S.&n-  0         so  ^JL  -  ■■ 


=  sX«o 


But,  it  is  possible  that  S  could  be 


moving  with  respect  to  the  fixed  frame  so  the  total  movement  is  given 

by: 

§  X  CU      4-  cfs 


Tela 


Ir 


If  there  exists  a  magnetic  field  in  the  direction  of  the 


Z  axis,  the  Larmor  precession  is: 


M 


it 


X(sxH 


Fig.  6.   Fixed  Reference  Frame 
Hence,  the  total  apparent  rotation  is: 


(  tfL  =  « » 


?-* 


But,  note  from  the  equations  concerning  the  Larmor  precession  that 
this  is  equivalent  to  a  magnetic  field  in  the  rotating  frame  which 

Heg  --  H  ■ 


has   the   value 
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Now  consider  a  pumping  light  modulated,  by  a  square  wave  at 
the  Larmor  frequency,  which  is  oriented  along  the  X  axis  of  a  fixed 
frame.   For  simplicity,  let  it  be  assumed  that  the  light  is  "on"  only 
when  the  X  axis  of  the  fixed  frame  is  aligned  with  the  XT  axis  of 
the  rotating  frame: 


2 


light 


%> 


? 


.^s 
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Ik 


observer 


Rotating 


Fixed 


Fig.  7.   The  Precessional  Reference  Frames  with 
Modulated  Light. 

If  the  angular  velocity  is  such  that  there  Is 

V 

no  residual  field  in  the  rotating  frame .   Further,  the  conditions  are 
the  same  if  the  observer  were  standing  on  the  rotating  frame  in  a 
darkened  room  with  the  pumping  light  blinking  in  synchronization  with 
the  rotation „   The  observer  sees  no  magnetic  field  and  no  movement 
of  the  spin  vector  in  the  rotating  frame  of  reference.   This  situation 
describes  a  static  condition  of  no  torque,  and  hence,  there  is  a 

minimum  of  light  absorption  which  might  be  ascribed  to  the  magnetic 

7 
field  effects. 

7cf,  J.R.  Singer,  (5),  pg.  48 
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A  closer  examination  of  this  resonant  condition  will  reveal  a 
space  quantization  along  the  X'  axis  which  is  due  to  the  pumping 

o 

light  alone.   This  effect,  first  suggested  by  Kastler  ,  was  later 
established  experimentally  by  Hawkins  (7)  using  sodium  atoms.   In 
effect,  the  space  quantization  is  the  same  as  the  magnetic  effects 
previously  described,  i.e.,  there  exists  preferred  orientations  of 
spin  vector  in  relation  to  the  X'  axis  in  the  presence  of  the  pumping 
light  and  in  the  absence  of  magnetic  field.   It  has  been  shown  in 
Appendices  A  and  B  that  transitions  involving  a  unit  change  in  angular 
momentum  are  associated  with  polarized  light.   If  the  spin  vector 
position  was  defined  to  be  positive  when  pointed  towards  the  observe  r, 
an  absorption  process  would  occur  when  the  vector  points  towards  the 
light  source.   Conversely,  when  the  vector  points  towards  the  observer, 
minimum  energy  is  extracted  from  the  beam.   Under  these  conditions, 
the  spin  vector  is  pointed  towards  the  observer  at  all  times  when 
resonance  is  established  in  the  rotating  frame.   This  is  the  equili- 
brium "pumped"  condition  of  the  sample,  and  no  transitions  are 
occurring  If  relaxation  effects  are  ignored. 

The  "off- resonance"  conditions  immediately  establish  a  residual 
magnetic  field,   Ht?<s?   oriented  parallel  to  the  Z  axis  which  has  a 

magnitude  He  —  UJ.      .   The  spin  vector  precesses  about  this  residual 
Y 

field,  bringing  it  into  position  to  be  "pumped"  back  to  its  equili- 
brium position.   Absorption  caused  by  the  pumping  process  is  due  to 
partial  mixing  among  the  Zeeman  sublevels,  and  if  fjg^   becomes 
very  large,  complete  mixing  will  occur. 

8A.  Kastler,  Proc .  Phys .  Soc .  A6_7,  853  (19  54). 
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It  can  be  assumed  that  the  spin  vector  motion  during  this 
woff- resonance"  condition  will  transcribe  a  complex  pattern  in 
space .   However,  a  simplified  picture  can  be  obtained  by  averaging 
random  instantaneous  "snapshots"  of  the  process  by  displaying  them 
in  a  continuum.   The  result  in  the  vicinity  of  resonance  might  show 
an  "expectation  cone"  as  the  surface  of  movement,  which  is  oriented 
towards  the  observer: 


4r-  L&rje     He% 


*^"sim3.tt"  fieg 


>  > £> 


A/o  Heg 


Fig.  8 „   Expectation  Position  of  the  Spin  Vector  in  a 
Rotating  Frame 

It  can  be  noted  that  the  net  effect  is  not  unlike  the 

new-standard  resonance  conditions  using  an  externally  applied  RF 

field  at  right  angles  to  the  static  field. 
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2.   The  Phe no mo logical  Approach 

Since  a  net  magnetic  moment  exists,  it  is  convenient  to  con- 
sider a  more  quantitative  viewpoint  of  the  resonance  phenomena. 
To  do  this,  the  individual  atomic  vector  parameters  may  be  com- 
bined to  give: 

M       -   The  net  instantaneous  resultant  magnetic  moment 
per  unit  volume  which  is  parallel  to  the  angular 
momentum  vector. 

A        -   The  instantaneous  resultant  angular  momentum  per 
unit  volume. 

7"   -  The  total  torque. 

f-j    -  The  static  magnetic  field  intensity. 

The  classical  relation  between  these  is: 

On  the  atomic  scale,  by  definition: 

where:   ^  -  gyromagnetic  ratio 

jU    ~  atomic  magnetic  moment 
^2  -  atomic  angular  momentum 
Hence,  by  analogy: 

Combining  the  above  relations  gives  a  result  which  is  similar  to 
the  previously  considered  precession  equations: 

^5  =  yj/ufij  a) 
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Let  the  magnetic  field  be  parallel  to  the  Z  axis.  The  major 
portion  of  the  precession  energy  is  due  to  the  magnetic  field  and 
can  be  written: 


t  =  -H 


o  >  *2 


Energy  changes  due  to  thermal  perturbations  will  be  then  caused 
by  changes  in  M%        -     Let  Hc         represent  the  equilibrium  value  of 
If  M^  d  j^c       ,   it  will  approach  this  value  exponentially  with  a  time 
constant  T;       which  is  the  "thermal"  or  "longitudinal"  relaxation 
time.   Therefore: 

Mb    _  -(tf*-M0) 

^  77    '  (2) 

Processes  which  do  not  effect  the  total  energy  can  be  attributed 
to  an  "effective  irregularity"  in  the  magnetic  field  of  order  fnf 
These  effects  can  be  characterized  by  a  "transversal  relaxation  time": 

T*  =  Tif' 

Assuming  that  the  transverse  effects  are  essentially 
exponential  in  character: 


and  that  these  effects  can  be  summed: 


(3) 
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}'■ 


Y[ft*H]..    +  **.  =o 


These  equations  define  the  motion  of  the  macroscopic 
magnetic  moment  with  the  magnetic  field  in  the  Z  direction, 
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3.   The  Optical  Modulation 

The  previous  phenomo Logical  equations  were  developed  under  the 
simplifying  conditions  that  the  equilibrium  position  coincided  with 
the  magnetic  field.   However,  from  the  intuitive  approach  presented 
previously,  it  is  obvious  that  the  equilibrium  position  for  the  pres- 
ent case  is  in  the  direction  of  the  X  axis  for  observation  along  the 
X  axis.   The  relaxation  factors  then  will  take  the  form: 


It  is  also  necessary  to  modify  the  equations  for  the  effect  of 
the  pumping  light.   Consider,  for  example,  the  effect  of  an  unpumped 
sample  being  irradiated  with  the  pumping  light.   The  rate  of  change 
of  one  sublevel  could  be  expressed  in  terms  of  the  number  of  atoms 
remaining  in  that  level  and  the  pumping  rate.   In  a  similar  manner, 
the  rate  of  change  of  the  X  component  of  the  net  magnetic  moment 
can  be  expressed  as  the  difference  between  the  final  magnitude  and 
the  instantaneous  magnitude  times  the  pumping  rate: 

(ilk)     =(hc-/-u)t\      (JjM      =  -/V*  (JmA 


Combining  these  relations  with  (1)  of  the  previous  section 

lfLMxHL   * 

5«a 


-  Y[m*hI    +  &  = 


n,  « 
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and  similarly  for  the  other  two  equations: 


o 


% 


(2) 


^^[flxffj.  +  £  -o  (3) 


where:  i  ( 

$>  '      b     ,  j_       $2   ~  p.   .+  j_ 

A         7J  £  >A     =   Pumping  rate 

These  are  the  equations  of  motion  of  the  net  magnetic  moment 
for  a  pumping  light  oriented,  at  right  angles  to  the  magnetic  field 
and  are  attributed  to  Bell  and  Bloom  (9)  . 

Since  the  precession  will  be  around  the  magnetic  field,  the 
net  XY  magnetic  moment  will  be  rotating  circularly  in  time.   It 
is  of  interest  to  consider  this  time  relation,  so  let: 

F  -  M*  +1  M^  |4) 

Inserting     ft  =  /fc    <X^  and  combining    (1)    through    (4)  : 

j£.    •+  I  *Hc  F    -r  £__   -,    /% 

It  s^  ^ 

Jr         T~  (6) 
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Due  to  the  light  modulation,  He   will  be  a  function  of  time; 
so  the  solution  of  (5)  must  be  restricted  by  defining^ 

-  The  polarization  that  would  be  produced  in  the  absence 
of  magnetic  field  and  thermal  relaxation  effects, 
such  that: 

Equation  (5)  is  then: 

g+*YH.r  +  £-*fl  (7) 

Now  let  the  light  be  modulated  by  defining; 

i  A 

j_  ^  ft  ^.  /  e#t  uyc 

Looking  for  solutions  of  the  form: 

and  separating  like  components  of  interest: 

L 
The  only  terms  which  give  an  observable  magnitude  are  those  associated 

with      since  at  resonance: 
Hence: 


\ 
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where:    £cO-Xlh'-tiJ 

From  the  previous  discussion  of  the  pumping  process,  it  can 
be  concluded  that  the  light  absorbed  would  be  proportional  to 
/J//-^7y-J     so  taking  the  real  part  of    : 


St         /f5        «* 


Using  a  sweep  field  frequency  which  is  slow  compared  to  the 
precession  frequency  and  a  detector  insensitive  to  the  precession 
frequency  will  result  in  a  signal  proportional  to: 


For  a  light  direction  vector  lying  at  the  angle      from 
the  pole  axis  (magnetic  field)  ,  the   ,'^   and  M »    above  will  con- 
tain a  sin  J-    dependence  resulting  in  the  signal  strength  correction 

p 
sin  9.  In  addition,  if  the  light  modulation  index  is  designated 

the  signal  will  take  the  form: 

5X  °L  — -- 

y  [/?*  +  c^f] 

It  can  be  seen  that  the  signal  will  increase,  with  the  square 
of  the  modulation  index;  the  light  intensity  and/or  the  sine  of  the 
angle  between  the  magnetic  field  and  the  light  directions . 
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4 .   The  Line  Shape 

For  a  given  set  of  operating  conditions,  the  signal  equation 
can  be  written: 


5  »    -^r  S  = 


A  /? 


If  the  frequency  is  varied  through  the  line  center  ( ¥ H0  ) 

the  one-half  amplitude  line  width  is  defined  by: 

r&  -  "u  =  ±  / 

and  therefore,  the  line  width  is: 
From  the  original  definition  of  /^    : 

.-t,     t-  where:   d2  -  The  average  pumping  rate 

P  *4 

ji    -   The  metastable  relaxation 
time 


7p    -   Average  pumping  time 


the  line  width  would  be: 


AUJL    -3L  (-L-    +  4r 
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These  times  are  difficult  to  predict  for  a  given  working 
situation.  Colgrove  and  Franken  (2)  have  measured  and 
found  it  to  be  about  one  quarter  of  a  millisecond  for  a  fixed 
condition  of  pumping  light  intensity,  absorption  cell  size  and 

pressure.   In  addition,  this  time  was  measured  using  unpolarized 

9 
light.   Phelps  and  Molnar  have  investigated  the  metastable  life- 
time in  some  detail  and  have  measured  lifetimes  as  high  as  20 
milliseconds.   From  this  work,  a  reasonable  estimate  for  the  pre- 
sent case  might  be  in  the  order  of  10"   seconds  when  the  cell 
pressure  is  0.5  mm  Hg.   If  it  can  be  assumed  that         is  of 
the  same  order  of  magnitude,  the  approximate  line  width  would  bes 


9 
A.Vo  Phelps  and  J.  P.  Molnar,  Lifetimes  of  Metastable  States 

of  Noble  Gases,  Phys .  Rev.  89,  1202,  (1953). 
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Ill .   The  Theoretical  Magnetometer  System 

A.  General  Design 

Instead,  of  sweeping  the  resonance  oscillator  to  obtain  a  signal 
as  in  previous  examples,  it  is  usually  more  convenient  to  modulate 
the  magnetic  field  at  a  slow  (audio)  rate.   The  comparison  of  the 
resulting  signal  with  the  sweep  field  cycle  then  allows  a  simple 
loop  control: 
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Fig.  9.   The  Theoretical  Magnetometer  System 
To  obtain  a  bright  pumping  light  source,  the  lamp  is  driven  by 
a  high  frequency  exciter  which  is  modulated  in  some  manner  at  the 
resonance  frequency.   The  exact  exciter  frequency  is  not  critical 
to  the  operation,  provided  it  can  be  modulated  at  1.^  mcs. 
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The  absorption  cell  is  also  provided  with  a  radio  frequency 
source  to  cause  ionization  with  subsequent  production  of  metastable 
atoms.   Again,  the  exact  frequency  is  a  matter  of  convenience,  but 
less  power  is  required  at  higher  frequencies. 

The  photodetector  can  be  any  type  sensitive  to  the  near  infra- 
red spectrum  which  can  respond  to  a  frequency  twice  that  of  the 
sweep  field.   The  lead  sulphide  detector  is  well  suited  to  this 
application  if  the  sweep  frequency  is  selected  below  about  1000  cps 

In  operation,  the  change  in  light  intensity  associated  with 

the  ambient  magnetic  field  causes  a  change  in  the  photodetector 

output  each  time  the  sweep  field  passes  through  zero.   This  signal 

is  amplified  and  compared  to  the  phase  of  the  sweep  field  to  detect 

a  system  error.   If  an  error  exists,  a  corrective  voltage  is 

developed  which  causes  a  corresponding  change  in  the  resonance 

oscillator  frequency  to  match  the  new  condition.   In  diagrammatic 

form: 

H0 -*-«s 


Field 


Signal 


Fig.  10.   Phase  Detector  Operation  when  Uf—KHo 
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When  the  instantaneous  frequency  does  not  match  the  ambient  field, 
the  signal  peaks  will  be  shifted  to  give  a  resultant  direct  current 
output: 


Field 


Signal 


tzfft  r 


Fig.  11 .   Phase  Detector  Operation  when  W  -f  JCHo 
This  is  an  over-simplified  representation  of  the  actual  opera- 
tion, since  it  is  assumed  that  the  signal  is  inserted  into  the  phase 
detector  as  a  polarized  pulse.   For  narrow  line  width  compared  to 
the  magnetic  sweep  magnitude,  A  OJjc   <<  Y  //5     this  is  not  a  serious 
error,  since  then  the  signal  will  appear  as  a  short  pulse  on  an 
otherwise  linear  sweep. 
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B.   The  Consideration  of  Noise 

Since  the  phase  detector  is  the  loop  error  control  device,  it 
is  instructive  to  examine  its  operation  in  some  detail. 

The  processing  of  the  raw  signal  could  conceivably  take  many 
forms  from  linear  amplification  to  actual  construction  of  a  new 
signal  from  some  correlation  device  inserted  between  the  photo- 
detector  and  the  phase  detector.   However,  the  phase  detector  opera- 
tion is  generally  quite  indifferent  to  the  average  shape  of  the  operat- 
ing signal,  and  the  experimental  system  complexity  should  be  kept  to  a 
minimum  consistant  with  adequate  performance .   In  light  of  these  con- 
siderations, the  signal  processing  must  consist  ,of  simple  amplification 

« 

with,  perhaps,  limiting  to  remove  variations  in  the  signal  amplitude. 
The  following  assumptions  form  a  model  quite  similar  to  the  experi- 
mental  conditions: 

a.  The  signal  is  limited  to  such  an  extent  that  a 
polarized  trapezoidal  pattern  results. 

b.  The  trapezoidal  sides  have  the  same  slope  as  the 
original  signal  at  the  one-half  amplitude  point. 

c.  An  "ideal"  phase  detector  is  used  which  has  the 
following  form: 
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Fig.  12.   Signal  Processing  Chain 


29 


Let  the  chopper  transient  occur  at  the  zero  crossover  of  the 
sweep  field  cycle  as  follows: 
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(No  Error) 


Fig.    13.      Phase   Detector  Waveforms 
If  the   current-to-voltage   convertor  is   a  capacitor  and  the 
signal  error  is   small  compared  to   the    trapezoidal  width: 


aE*  =  tie  ATe 


(i) 


where:  & £c   is  the  error  voltage  output  per  sweep  field 
cycle . 

£le   is  the  time  displacement  of  the  processed 
signal  from  the  mean. 

Due  to  the  severe  limiting  of  the  signal,  the  input  noise 

fluctuations  will  only  result  in  a  displacement  of  the  trapezoid: 

30 


-Li 


Amplitude  limit 


a*   TV 

Fig.  14-.   Noise  Variations  in  the  Processed  Signal 
where:  l7~      is  the  rms  value  of  the  input  noise 

&t     is  the  rms  value  of  the  noise  displacement 
TV   is  the  time  of  rise  of  the  trapezoid 
The  slope  of  the  unprocessed  signal  at  the  one-half  amplitude 
point  can  be  found  by  differentiation  of  the  signal  equation  and 
substitution  of  the  relation: 


(  ^    XJk  -J*  =  / 


Hence: 


s  - 
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/M 


**+  a /v* -<-)"• 
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But  the  one-half  amplitude  line  width  was  previously  determined 
to  be  £  (^   ,  so 

J£ 

slope^  "         AUJJL 
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The  field  is  being  swept  at  a  low  audio  rate.   If  this 
corresponds  to  a  resonance  frequency  rate  of  vJfo    radians/sec  , 
the  slope  of  the  signal  in  time  will  be: 

slooe  =   —  W 
and  the  relation  with  the  processed  signal  is: 

ASK    _     BKWa. 

where  the  dimension  change  $  K   amperes  is  necessary  to  maintain 
a  constant  amplitude  ratio  and  like  units.   From  the  similar 
triangle  relations  of  the  trapezoid: 
At  ~  1*   -  ^uJt 


3  tr-       XP  3K  Wk 


The  rms uncertainty  in  signal  position  due  to  noise  is 
therefore: 


From    (1) 


Ate    '    ±AiI!i 

c  Wh 

And  the  "rms  error  voltage  per  unit  signal  per  cycle"  is: 

Although  further  processing  will  modify  the  actual  error  due 
to  noise,  the  line  width/signal-to-noise  relation  will  determine 
the  input  error.   The  conclusion  is:   The  magnitude  of  the  noise 
error  is  proportional  to  the  line  width  and  inversely  proportional 
to  the  signal-to-noise  ratio. 


32 


IV.   The  Experimental  Magnetometer  System 
A.   Principles  of  Construction 

The  guiding  principles  for  construction  of  the  magnetometer  were: 

1)  Obtain  a  low  line  width/signal- to~noise  ratio  consistant 
with  local  conditions. 

2)  Use  standard  equipment  where  possible. 

3)  Provide  variable  filtering  and  system  gains 
to  optimize  the  field  tracking  ability. 

To  optimize  the  line  width/  signal-to-noise  ratio  there  must  be 
homogeneity  of  the  magnetic  field  within  the  absorption  cell  and  an 
orientation  of  the  optical  system  perpendicular  to  the  local  field. 
The  magnetic  field  direction  near  the  system  location  was  found  to 
be  33  degrees  from  the  vertical.   Previous  experiments  had  indicated 
that  a  maximum  signal-to-noise  ratio  of  approximately  ten  could  be 
achieved  in  the  vertical  position,  and  the  optical  system  was  so 
mounted,  although  this  represented  a  theoretical  signal  amplitude 
which  was  30%  of  that  obtainable  in  the  optimum  orientation.   This 
decision  considerably  eased  the  optical-structure  details.   The 
field  homogeneity  was  dictated  by  local  working  space,  the  equipment 
location  and  orientation,  and  the  total  field-shaping  efforts.   Efforts 
were  made  to  use  non-magnetic  components  in  the  immediate  vicinity  of 
the  cell  (a  few  feet)  ,  but  beyond  this  no  particular  efforts  were 
made  to  improve  the  gross  field. 
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B.   The  System  Description 

Following  the  general  design  of  the  theoretical  magnetometer, 
the  experimental  equipment  was  assembled  into  the  form  shown  in 
the  following  block  diagram: 
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Fig.    15.      The   Experimental  Magnetometer  System 
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Due  to  the  high  stray  60  cps  magnetic  fields  found  in  the 
area,,  no  additional  sweep  field  was  used.   The  magnitude  and  phase 
of  this  field  could  be  adjusted  slightly  by  proper  orientation  of 
major  components  using  60  cps  power-primarily  power  supplies. 

The  major  commercial  equipments  used  in  the  system  were: 

V.A.   354-  kc.  Oscillator 

V.A.   Mixer  and  0-3  kc  Filter 

V.A.   G-10  Recorder 

H.P.   Frequency  Meter  500B 

H.P.   Amplifier  450A 

Krohn-Hite  Filter  330A 

Kodak   Ektron  Photoconductive  Detector 

The  schematic  drawings  of  the  remaining  equipments  are  in- 
cluded in  Appendix  D,  with  the  exception  of  the  30  mc  exciter  and 
the  isolation  amplifier. 

Filament  power  for  the  modulator,  the  115  mc  exciter  and  the 
30  mc .  exciter  was  provided  by  battery,  and  the  high  voltage  supply 
was  provided  from  commercial  power  supplies.   The  preamplifier, 
frequency  control  unit,  1 . M-  mc.  oscillator,  harmonic  generator  and 
recorder  bias  were  completely  battery  supplied. 

Photographs  of  the  working  system  are  shown  in  Figs.  16,  17 
and  18. 

The  heart  of  the  magnetometer  system  is  the  aluminum  optical 
rack  shown  attached  to  the  wall.   As  can  be  noted  from  the  interior 
photograph,  the  helium  lamp  and  the  115  mc .  exciter  occupy  the  low- 
est tier  together  with  a  spherical  mirror  which  increases  and 
collimates  the  usable  pumping  light.   The  light  then  passes  through 
a  linear  polarizer  (Polaroid  HR)  and  a  quarter  wave  plate,  forming 
circularly  polarized  infra-red  light  components  before  entering  the 
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absorption  cell.   The  absorption  cell  is  excited  into  a  weak  dis- 
charge by  a  30  mc .  oscillator  connected  to  a  pair  of  wire  rings 
wrapped  around  the  cell.   (It  must  be  noted  at  this  point  that 
the  brilliance  of  the  cell  was  increased,  and  that  of  the  lamp 
decreased  from  optimum  in  order  to  obtain  a  photograph.)   The 
light  not  absorbed  in  the  cell  continues  through  a  short  focal- 
length  lens  system  (not  shown) ,  an  infra-red  filter  to  reduce  ex- 
terior interference  and  then  strikes  the  photodetector  contained 
in  the  preamplifier  unit.   As  the  stray  60  cps  magnetic  field 
changes  the  total  field  within  the  cell,  a  decrease  in  the  light 
absorption  occurs  when  U)0   -  )(  He  •   This  is  registered  as  a 
change  in  resistance  at  the  photodetector . 

The  remainder  of  the  equipment  processes  the  signal  and 
develops  a  direct-current  error  signal  which  is  applied  to  a 
voltage-variable  capacitor  to  control  the  1.4-  mc .  oscillator  for 
frequency-locking  on  the  line  center.   The  system  gains  are  shown 
on  the  block  diagram. 

Modulation  of  the  lamp  is  accomplished  by  plate  and  screen 
modulation  of  the  115  mc .  exciter.   An  improvement  in  the  signal 
was  noted  when  the  modulation  of  the  exciter  exceeded  100%, 
probably  indicating  that  the  actual  light  modulation  is  less 
than  that  of  the  exciter. 

The  magnitude  of  the  field  changes  can  be  converted  directly 
from  the  oscillator  frequency  by  application  of  the  factor  28 
cycles/  gamma.   To  record  the  frequency  variations,  the  fourth 
harmonic  of  a  commercial  crystal  oscillator  was  mixed  with  the 
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system  1.4  mc  oscillator  signal  to  form  a  difference  frequency 
in  the  audio  range.   This  in  turn  is  converted  to  a  direct-curr. 
signal  which  is  proportional  to  the  difference  frequency  and 
applied  through  a  calibration  network  to  a  pen  recorder. 

The  major  components  of  the  system  are  noted  on  the  photo- 
graphs, including  the  field- shapers  which  are  used  to  modify  the 
field  homogeneity  within  the  cell  or  to  adjust  the  total  field  to 
lie  within  the  audio  range  for  recording  purposes. 
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The  Experimental  ha^ne to meter 
-  Left  Hand  View 

1.  Optical  Rack 

2.  Detector  and  Preamplifier 

3.  Absorption  Cell   30  inc.   Exciter 

4.  1.4  mc   Modulator 
b.    1.4  mc.   Oscillator  and 

Frequency  Control   Unit. 
6.   Field-Shaper 
'7.    20  db  Amplifier 

8.  Active  Filter 

9.  Phase  Detector 

10.  harmonic  Generator 

11.  Power  Supplies 
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Ti.e  Exrerimental  Magnetometer 
-  Ri&ht  hand  View 
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V.   The  Experimental  System  Performance 

A.   The  Measurements  of  the  Magnetic  Field 

The  helium  magnetometer  system  was  designed  primarily  to 
measure  changes  in  the  magnetic  field.   Theoretically,  the  total 
absolute  magnetic  field  can  be  obtained  from  the  sum  of  the  crystal 
frequency  and  the  difference  frequency  produced  by  mixing  the  reson- 
ance oscillator  and  crystal  frequencies.   However,  in  the  present 
case,  only  a  portion  of  the  difference  frequency  is  used  by  the  re- 
corder.  The  portion  of  the  difference  frequency  used  is  selected 
by  the  recorder  biasing  and  calibration  network.   To  further  reduce 
noise  perturbations  in  the  recording  system,  a  low  pass  filter,  which 
has  a  3  kc  upper  cut-off,  restricts  the  usable  difference  frequency 
to  a  range  3  kc  either  side  of  the  crystal  frequency.   Two  factors, 
therefore,  are  apparent.   1)  The  total- field  information  is  lost  in 
the  recording  system,  and  only  magnetic  field  perturbations  will  be 
recorded.   2)   The  magnetometer  precession  frequency,  as  taken  from 
the  resonance  oscillator,  must  be  within  3  kc  of  the  crystal  standard 
for  recording  purposes. 

To  record  magnetic  field  perturbations,  the  magnitude  of  the 
total  field  within  the  cell  was  controlled  by  use  of  field-shaping 
metal  objects  which  were  placed  near  the  helium  system.   These 
objects  generally  cause  magnetic  field  inhomogeneity  within  the  cell 
which  increases  the  helium  line  width.   As  previously  shown  under 
the  consideration  of  noise,  an  increase  in  the  line  width  causes  a 
degradation  of  the  performance  due  to  the  contribution  of  noise.   In 
the  present  case,  this  was  considered  acceptable  in  order  to  obtain 
recordings  of  the  magnetic  field. 
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The  general  performance  of  the  experimental  device  was 
measured  by  simultaneous  recordings  with  a  commercial  magneto- 
meter which  used  rubidium  as  the  sensing  element.   The  rubidium 
sensing  head  was  located  in  an  open  field  approximately  200  feet 
from  the  helium  magnetometer  location. 

The  first  closed-loop  operation  of  the  magnetometer  is 
shown  in  Fig.  19.   This  recording,  taken  between  1900-2400  GMT 
on  27  February,  1961,  shows  the  local  mid-day  diurnal  change  as 
recorded  by  the  rubidium  (left)  and  the  helium  (right) .   The 
time  scale  increases  from  bottom  to  top.   The  helium  system  was 
operating  on  one  side  of  its  crystal  standard  and  the  rubidium 
system  on  the  opposite  side  of  its  crystal  standard;  so  ideally., 
the  two  traces  should  be  mirror  images .   Due  to  the  experimental 
nature  of  the  helium  system,  a  certain  degree  of  frequency  drift 
was  experienced  in  this  recording.   The  broadening  of  the  helium 
trace  is  attributed  to  damped  oscillations  of  the  servo-loop.   The 
primary  objective  of  this  recording  was  to  test  the  system  ability 
to  track  a  slow  and  relatively  predictable  change  in  the  magnetic 
field.   A  smaller  change  was  also  recorded  simultaneously  by  both 
systems  (see  Fig.  19)  . 

The  second  recording  is  shown  in  Fig.  20.   The  upper  trace 
pair  is  a  time-continuation  of  the  lower  trace  pair,  and  the  helium 
system  output  is  the  upper  trace.   Disturbances  between  1500-1530 
were  caused  by  movement  of  the  operator  within  the  room  containing 
the  magnetometer.   The  disturbance  near  1730  was  caused  by  opera- 
tion of  electrical  equipment  in  a  nearby  building.   An  approximate 
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two  gamma  shift  was  experienced  at  17  20  due  to  changes  in  the 
magnetic  sweep  (60cps) ,  so  the  graph  was  broken  at  this  point 
and  moved  upward  to  match  the  general  field  movement.   Signal- 
to-noise  decrease,  due  to  decaying  batteries  within  the  reson- 
ance oscillator,  progressively  degrades  the  helium  performance. 
This  is  particularly  noticable  between  2000  and  2100  when  the 
pumping  light  intensity  was  increased  to  improve  the  signal-to- 
noise  ratio.   Note  that  the  quality  of  system  operation  approaches 
that  of  the  rubidium  magnetometer  in  the  period  1610-1700,  al- 
though the  minor  perturbations  within  this  period  may  not  have 
affected  each  magnetometer.   Disregarding  the  man-made  disturb- 
ances, the  helium  tracking  was  within  _+  3  gamma  of  the  commercial 
magnetometer  over  the  M-5  gamma  change  which  required  4-  hours  to 
record. 

The  third  recording,  Fig.  21,  shows  the  relative  system 
tracking  ability  for  major  field  changes  which  occurred  over 
short  periods  of  time.   The  difference  in  scales  is  attributed 
to  the  recording  equipment.   In  general,  the  helium  system  was 
within  0  5  gamma  of  the  commercial  magnetometer  over  short  periods 
of  time. 

The  fourth  recording  pair,  Figs.  22  and  23,  show  a  large 
perturbation  in  the  field.   Again  the  scale  discrepancy  is  in 
the  recording  equipment.   It  is  noteworthy  that  many  minor  field 
perturbations  of  the  order  of  a  few  tenths  of  a  gamma  were 
recorded  simultaneously  by  the  magnetometers. 
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From  the  graphs  obtained,  it  can  be  concluded  that  the 
helium  magnetometer  was  operating  within  0.5  gamma  of  the  com- 
mercial magnetometer  over  short  periods  of  time,  (several  minutes) 
and  under  optimum  conditions.   Due  to  drifting  characteristics 
caused  by  battery  decay,  no  approximation  can  be  made  for  long- 
term  accuracy,  although  the  diurnal  correlations  prove  tracking 
ability. 

Unfortunately,  the  term  "optimum",  as  used  in  this  section, 
means  "optimum  considering  the  circumstances  and  environment"; 
and  therefore,  the  true  accuracy  under  conditions  similar  to  the 
commercial  magnetometer  environment  cannot  be  determined  from 
the  data  obtained. 
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Fig.  19  System  performance-  27  February,  1961 
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Fig.    20  System  performance-   4  March,    1961 
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Fi&.  21  System  performance-  p-6  March  I96I 
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Fig.  22  Large  Field  Perturbation-Rubidium  Magnetometer 
o  March,  I96I 
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Fig.  23  Large  Field  perturbation-  Helium  Magnetometer 
6  March,  I96I 
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B.   The  Subsidiary  Measurements 
1„   The  Absorption  Phenomena 

Although  the  analysis  of  the  system  operation  is  not  depend- 
ent on  a  knowledge  of  the  magnitude  of  the  absorption,  it  is  a 
graphic  verification  of  the  theory.   To  obtain  information  con- 
cerning the  magnitude  of  absorption,  photographs  were  taken  of 
the  photodetector  voltage  variations  when  the  absorption  cell 
discharge  was  turned  on  and  off  in  the  presence  of  and  in  the 
absence  of  the  pumping  light.   The  following  equipment  config- 
uration was  used: 
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Fig.  2^.   Equipment  Configuration  for  Absorption  Observation 
In  the  figures  to  follow,  a  light  intensity  increase  is  re- 
corded as  an  increase  in  the  voltage  traces  (i.e.  upward  deflection 
of  the  trace)  and  slow-changing  transients  are  attributed  to  RC 
charging  times  within  the  equipment. 
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In  Fig.  25  an  increase  in  the  light  intensity  at  the  detector 
was  recorded  when  the  pumping  light  remained  on  and  the  cell  dis- 
charge was  turned  off.   The  magnitude  of  the  change  is  approximately 
600  millivolts  with  a  detector  bias  of  38  volts.   This  is  a  1.6% 
increase . 

Fig.  26  shows  a  decrease  in  the  light  intensity  at  the  de- 
tector when  the  pumping  light  remains  off  and  the  cell  discharge  is 
turned  off.   The  magnitude  ot  this  cnange  is  approximately  50 
millivolts  or  a  0.1%  decrease. 

The  signal  under  these  conditions  is  shown  in  Fig.  27.   The 
measured  signal-to-noise  ratio  was  four,  and  the  peak  amplitude 
was  approximately  600  uv .   The  resulting  signal- to-light 
absorbed  is  approximately: 

£f*>  tit,-f   ^  /X  to'3 
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Fig.  25  Detector  Voltage  Variation- Lamp  on,  Cell  off. 
Scale:  0.1  volt/div,  10  msec/div 
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Fig.  26  Detector  Voltage  Variation-Lamp  off,  Cell  off. 
Scale;  0.02  volt/div,  10  msec/div. 


Fig.  27  The  Unprocessed  Signal.  Scale:  50  mv/div 
2  msec/div. 
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2.  The  Signal  Processing 

The  average  signal  obtained  in  the  33  degree  orientation 
is  shown  in  Fig.  28.  The  signal-to-noise  ratio  was  measured 
to  be  seven.  By  filtering  to  a  100-240  cycle  bandwidth,  a 
considerable  improvement  could  be  achieved  aa  shown  in  Fig. 
29  as  a  Lissajous  and  Fig.  30  on  a  linear  sweep. 
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Fig.  28  Unprocessed  Signal      Fig.  29  Filtered  Signal 

S/N  7  Bandwidth  100-240  cps 


Fig.  30  Processed  Signal 

Bandwidth  100-240  cps 
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3 .   The  Measurement  of  Line  Width 

The  signal  entering  the  phase  detector  is  symmetric  in  time 
only  when  the  precession  frequency  coincides  with  the  resonance 
oscillator  frequency.   This  frequency  corresponds  to  the  ambient 
magnetic  field.   As  the  resonance  oscillator  is  changed  from 
this  "correct"  frequency,  the  phase  detector  will  develop  an 
error  voltage  which  has  a  polarity  determined  by  the  direction 
of  frequency  change.   The  magnitude  of  the  error  voltage  will 
continue  to  increase  until  the  frequency  error  reaches  one-half 
the  line  width,  at  which  point  the  servo-system  locking  will 
become  unstable.   By  use  of  a  frequency  counter,  these  points 
were  measured  under  several  representative  sweep  field  conditions 
The  results  were  averaged  to  give  an  order-of-magnitude  line 
width  of  approximately  10,000  cps. 
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VI .   Conclusions 

A.   System  Performance 

The  objective  of  this  work  has  been  to  show  that  a  modulated, 
light  helium  magnetometer  could  be  constructed  and  could  be  caused 
to  follow  the  changes  in  the  earth's  magnetic  field.   During  this 
endeavor,  no  particular  value  was  placed  on  accuracy,  per  se, 
provided  the  resulting  performance  could  be  clearly  associated 
with  changes  in  the  ambient  field.   The  actual  performance  con- 
firmed the  fact  that  the  system  met  this  specification.   During 
short  periods  of  time,  the  performance  of  the  system  was  within 
0.5  gamma  of  the  rubidium  magne tome te r „ 

Several  modifications  could  be  made  to  improve  the  accuracy 
of  the  experimental  model; 

1.  Orientation: 

The  orientation  used  was  such  that  only  30%  of  the 
theoretical  maximum  signal  amplitude  was  obtained.   By  obtaining 
the  maximum  signal  amplitude,  a  three-fold  increase  in  signal-to- 
noise  ratio  could  be  expected. 

2 .  Environment 

The  experimental  equipment  was  installed  in  a  building 
which  contained  a  variable  60  cycle  sweep  field  and  equipment 
which  caused  magnetic  interference.   These  factors  were  noted  to 
cause  as  much  as  two  gamma  shifts  in  the  recorded  output  together 
with  additional  noise  perturbations.   A  magnetically-clean 
environment  could  eliminate  these  problems. 
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3 .  Packaging 

The  system  components  used  were  bulky  and,  for  the  most 
part,  inadequately  shielded.   This  added  to  the  "magnetic- 
clutter"  previously  mentioned  under  environment,  and  it  also 
provided  numerous  sources  of  noise,,   By  designing  equipment 
specifically  for  the  system,  a  general  improvement  in  performance 
could  be  expected. 

4.  Line  Width 

The  line  width  directly  effects  the  accuracy  of  the  system. 
Magnetic  field  homogeneity  in  the  absorption  cell  will  determine 
the  line  width  for  a  given  operating  condition.   The  line  width 
measured  in  the  present  case  was  approximately  10,000  cps,  which 
is  considerably  more  than  one  might  estimate  from  a  theoretical 
value.   An  effort  to  clean  up  the  magnetic  environment  and  the 
addition  of  a  homogeneous  magnetic  sweep  would  result  in  greater 
system  accuracy. 

5.  Reference  Crystal  Oscillator 

The  recording  system  accuracy  is  dependent  on  the  accuracy 
of  the  reference  crystal  oscillator,  the  frequency  meter  and  the 
recording  equipment.   In  the  present  case,  a  crystal  reference 
oscillator  was  used  which  contained  no  shielding  from  other  equip- 
ment and  no  temperature  stabilization .   By  using  a  reference  system 
which  has  a  short  term  stability  of  one  part  in  a  million,  a  .05 
gamma  error  would  be  present  which  could  not  be  improved  by  changes 
in  the  magnetometer.   The  inaccuracies  of  the  frequency  meter  and 
the  recording  equipment  further  degrade  the  recorded  performance. 
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With  the  above  modifications,  the  experimental  performance 
could  be  expected  to  improve.   With  sufficient  care,  the  accuracy 
of  the  system  could  conceivably  approach  the  limitations  imposed 
by  the  reference  oscillator. 
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APPENDIX  A 
ZEEMAN  LEVEL  ENERGY  CONSIDERATIONS 
From  classical  mechanics,  the  force  exerted  on  a  magnetic 
monopole  is  given  by: 

Where:  .-w    =   magnetic  charge 

H     =  magnetic   field  intensity 
The    torque    acting  on   a  pivoted  monopole  will  be: 
T  a  I^JUn  0 


Where:     £  ~  length  of   the  pivot 


arm, 


Q   =  angle  between  the  monopole 
axis  and  the  field. 

But,^/^  is  simply  the  magnetic  moment,  ^U    so: 
T   =  M  H  -4**  & 

The  Energy  associated  with  this  torque  is: 
O 

*/  (1) 

where  £  =  0   corresponds  to  the  field- free  condition. 

Let  L  be  the  net  angular  momentum.   Then,  for  a  charge,  <** 
moving  in  a  circular  orbit  at  angular  rate  &>     : 

L   =  tv  r*  oJ 
But,         . 

^  =  ^  A 


£^A 


=  tA 

y  where  a,    =   the  filamentary  current 

established  by  the  rotating 
charge „ 

/  =  period  of  orbit.  =  "T~ 
A    -   area  of  orbit. 
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Hence: 


Expressing  the  charge  in  electrostatic  units  and  combining 
these  results  for  an  electron: 


The  angular  orbital  momentum  of  the  electron  is: 


(2) 


A  /r?Tw) 


3l  ir  v  -  ■  -  (3) 

snd  that  of  the  spin  has  been  found  to  be  twice  that  of 
(3) 


Combining  (3)  and  (M-)  into  (2)  : 


ML    - 


Ma 


-e 


4«r^ 


^cV/l(L+<)    *l    ='^/UL+'j^    ( 


5) 


2elltCJs(s+\)  2s     z-jMffJsCs+i)  a. 


(6) 


where  ,&&    =  Bohr  magneton  =  0.9273  x  10"  °  erg/oersted 
Drawing  these  results  in  vector  form: 


Angular 
momentum 


Magnetic 
moment 


Fig .31. Vector  sum  of  the  Moments  and  Momenta 
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It  can  be  noted  that  where  spin  is  involved  the  net 
angular  momentum  vector  is  not  aligned  with  the  net  magnetic 
moment.   The  component  of  the  net  magnetic  moment  which  is  then 
of  interest  can  be  calculated: 


s-tj-^L  <**t#*jJol  +s*s  **D**,>*/J 


(7) 


But  from  the  law  of  cosines 


2  lf(S+i)L(L  +  i) 


(8) 


and  combining  (5)  through  (9) 


=  {i *  yy;^aflj}/^)  -  ^  i/^>) 


V 
Aer        axcx-Hj  (io) 


where  ^   is  the  Lande  g- factor. 
Combining  (1)  and  (10)  : 

£  =  <f/ts  Jrcr+i)  \\c<x>e  (ii) 

which  can  be  expressed  in  terms  of  the  yj"~(<jr-w)    projection 
on  the  magnetic  axis: 

4E=  5.^h(am)  d2) 


where:  /^l  ar  jT(J~t-i)    UX)  & 
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APPENDIX  B 
THE  EFFECT  OF  POLARIZATION  ON  TRANSITIONS 
The  intensity  of  radiation  from  the  atom  is  dependent  on  the 
magnitude  of  the  change  in  the  dipole  moment.   In  the  case  of  two 
charges   £   and  —  £   positioned  a  distance  o      from  each  other, 
the  dipole  moment  is  a  vector  of  magnitude  So        •   In  an  ensemble 
of  point  charges,  the  dipole  moment  vector  could  be  defined  by  the 
cartesian  components; 

£-**xi        X€^       2><^ 

<<  L  L 

In  the  case  of  a  single  point  charge  whose  exact  position  is 
unknown,  a  certain  fraction  of  the  total  charge  £   can  reasonably 
be  considered  to  be  at  a  fixed  point  if  the  distribution  probability 
is  known.   Since  the  probability  density  function  is  Y   T  , 

from  quantum  mechanics,  the  electric  density  would  be: 

and  the  electric  dipole  moment  components  for  a  single  electron 
would  be : 


*  '<v 


=     €   %,    t,    X    J1 


Pv  =J*  r„  tirJr 


^Je  L  tfzJr 


(i) 


(2) 


(3) 
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where   J    is  the  element  of  volume  and  the  subscript  represents 
the  stationary  state  /?} 

If  the  transition  from  state  sn  to  state  swi  is  made, 
the  electric  moment  (transition  moment)  P  can  be  described  by 
the  same  equations  with  the  change    jt.-/i     I  for  I     \m 

It  may  be  recalled  that  the  Schrodinger  equation  can  be 
written: 

tl  +  ? T  +  til  +  v T  =  jj>    ,- T 
df        "37*        ^r*J 


II 


which  has  a  solution: 


f=  fe-J**®* 


CO 

where:  j    is  the  time  independent  solution, 

t      =   Energy 

£   =  time 

h  =  Planck's  constant 

Considering  the  time  factor  solution  applied  to  the  two  states, 

results  in: 

r~'22  It 


ajyfh  I  fc-.'W   fc^  1 


e 

(5) 

where  there  exists  the  familiar  relation  that  the  energy  difference 
is  related  to  the  frequency  by: 

h 

If    m      is  a  vector  which  defines  the  amplitude  of  the 

vibration  of  the  electric  moment,    P      then  the  combination  of 
(1)  ,  (2)  ,  (3)  and  (5)  results  in: 
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as  the  time-variable  electric  moment  associated  with  the 
transition  ./yi      to  .^w  .   f\  has  the  components: 


ft"-  ft,  UaJr  j    eve 


(6) 

Considering  a  magnetic  field  direction  along  the  Z  axis, 
the  characteristic  solutions  to  the  Schrodinger  equation  would 
take  the  form: 

(7) 
since  the  states  are  rotationally  symmetric  about  the  magnetic 
field. 

Combining  (6)  and  (7) : 

fC  w  =  fffs*eJ  ^  ~**s)  * *  ^ 

and  similarly  for  y       and   j7 

In  cylindrical  coordinates: 

jr~  CJcpJCcU 

y  —  ?    A^n   <p 
The  Z  component  of  \\  is  then: 

t£~=  fffj.3  *  eJ  M*  z  e  J*  J?  Jz 
=  ff*3*z e  Jz  Jtje*M*J* 

where:      Asm  =  /*>?£' 
This    integral   is   non-zero   only  when: 

**„  J  (**»)*  u 

6*  t  is  non-zero. 

An  examination  of  this  latter  integral  will  reveal  that  the  result 

will  not  be  zero  if;  /^.fTl   -r  O 
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This,  of  course,  indicates  that  only  those  transitions  involving 
=  0  can  occur  parallel  to  the  magnetic  field. 
The  x  component  of   |      is: 

R™=  fff6  y>Jz  Jc  feJ^'^tJ* 

where    the   latter  integral  can  be  written: 

■J-  (\e J (A" * ° *  +  e J (*"" ' '^ }J 4 

This  has  a  non-zero  value  only  if   <&>»1  =  4_  /   .   Similarly,  for 
the  y  component: 

_L_  f  reJ(*»*-»)<*  _ed  (*"-')  <*\  J4 

Now  let  /X/YA  =-1  .   Examining  the  4    dependence  of  the  x 
component: 

The  y  component  is: 

tf~  =  -j  c-(r,  z)  tr- 

Hence,    the   vector      /\  is: 

j?"""  =  F(t 2)  M  a*  -J  c-(C  z)  t  S^ 

Looking  into  the  Z  axis,  the  y  component  lags  behind  the  x 
component  by  90  degrees  in  time,  and  the  rotation  is  counterclock- 
wise.  An  investigation  of   A^tt  =  -f-  /   will  reveal  a  circular 
polarization  of  opposite  direction. 
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APPENDIX  C 
THE  POPULATION  DENSITY  CONSIDERATIONS 

As  described  previously,  the  general  pumping  action  is  an 
absorption  process  which  effectively  excites  the  individual  atom 
from  a  35  sublevel  to  a  ^P  sublevel.   Each  allowed  transition  has 
a  finite  relative  probability  of  occurrence  which  can  be  calcu- 
lated  ,  but  the  need  does  not  arise  here.   The  allowed  downward 
transitions  can  be  assumed  to  be  primarily  due  to  the  spontaneous 
emission. 

In  the  case  of  gas  pressures  in  the  order  of  a  few  millimeters 
or  higher,  there  is  a  considers  Jle  disorientation  and  loss  of 
phase  memory  in  the   P  levels  before  the  downward  transitions  can 
occur.   This  allows  the  simplifying  assumption,  often  used,  that 
all  the  -*S  sublevels  have  an  equal  probability  of  being  repopulated 
by  the  returning  atoms . (2)  (9)   In  the  present  case,  however,  the 
gas  pressure  is  in  the  order  of  one-half  millimeter,  and  it  must 
be  assumed  that  no  such  disorientation  exists. 

The  primary  interest  here  is  to  show  that  an  orientation  can 
exist  by  virtue  of  a  net  population  increase  of  one   S  sublevel 
at  the  expense  of  the  other  two  in  the  absence  of  thermal  relax- 
ation effects. 

The  allowed  helium  transitions  which  are  of  interest  are 
shown  in  Fig.  32. 

^Condon  and  Shortley,  The  Theory  of  Atomic  Spectra, 
Cambridge  Press,  193  5,  Chap.  3. 
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Fig.  32.   The  Allowed  Transitionsin  Helium 


A  table  of  upward-going  transition  rates  can  be  formed: 
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*      ifttfl  '  *.'   !  ^ 
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Table  1.   Absorption  Transition  Rates 

Only  the  circled  values  will  be  allowed  with  circularly  polarized 
light  in  accordance  with  Appendix  B,  where   AM  —  +  f  „ 

The  relative  spontaneous  emission  probabilities  have  been 
computed  by  Colgrove  and  Franken  (2)  : 
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Table  2 .   Relative  Spontaneous  Emission  Probabilities 

Comparing  these  two  tables,  a  tabulation  can  be  made  to 
replace  Table  2  when  using  polarized,  light: 
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Table  3„   Spontaneous  Emission  Probabilities  in 
terms  of  Absorption  Rates 


The  rate  equations  which  describe  the  change  in  population 

densities  may  now  be  written.   Let  C  K+  gi-ve  the  absolute 

1 

absorption  rate  for  the  ^s,  +1  level,  etc.   Then: 


c    IT  " 

J...  a  Wo 


CW/*  —   -.  .—  ^ 


\f  ?el  + 


*■  *-j 
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where  sff.        ^   and  ,<??     denote  the  instantaneous  sublevel 
populations  and   /U  ■=:  sri  +  -f  ^c  *-/n~    • 

Note  that  the  atoms  pumped  from  the  +1  sublevel  return  to  that 
level,  whereas  the  other  two  levels  are  being  depleted.   Assuming 
equal  pumping  light  intensities  and  no  relaxation  effects,  the  +1 
sublevel  population  will  become   ^. 
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APPENDIX  D 
THE  EQUIPMENT  SCHEMATIC  DIAGRAMS 
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Diagram  1  Photodetector  preamplifier 

71 


52 


Diagram  2  1.4  mc,  Oscillator  and  Frequency  Control  Unit 
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Diagram  3     1.4  mc.    Modulator  and  115  mc.   Lamp  Exciter 
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Diagram  4  60  cycle  Phase  Detector 
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